天然气水合物体系的相平衡及甲烷溶解度计算          下载免费PDF全文

苏正  陈多福  冯东  刘芊 《地球物理学进展》2007,22(1):239-246

海洋环境中天然气水合物的形成除了合适的温压条件外，还必须有充分的甲烷供给.本文介绍了甲烷-水体系的甲烷饱和溶解度、水合物体系中甲烷水合物溶解度计算方法.在气-液二相平衡甲烷饱和溶解度计算中，关键在于状态方程的选择和合适的混合规则的运用，Duan的计算模型在温度、压力和盐度变化上都具有很大的适用性，且易于应用.在含水合物的相平衡体系中，在已知组分和假定可能存在相的前提下，可利用模拟退火算法优化总吉布斯自由能，确定是二相还是三相体系，并求解甲烷水合物溶解度.在海水环境下盐的存在使平衡发生移动，利用德拜—休克尔理论或Pitzer电解质溶液理论校正盐度对于海水活度的影响，求解海水环境中甲烷水合物溶解度.基于气-液二相平衡理论的K-K方程，在临近水合物生成条件下实验或计算确定亨利常数等参数后，可计算三相平衡甲烷水合物溶解度，且简单易用.  相似文献   

伊利石结晶度Kübler，Weaver和Weber指数关系式常数的确定与应用   总被引：1，自引：0，他引：1  

王河锦  陶晓风Rahn Meinert 《高校地质学报》2007,13(3):561-565

根据衍射峰具有的不同分布特点推导出了14个伊利石结晶度Kübler指数，Weaver指数和Weber指数关系式中常数项C1-C14的数值，提出了PearsonⅦ分布中形态参数μ与C5，C6，C11，C12的4个定量方程。由此，14个伊利石结晶度指数关系式可直接用于实际计算。根据阿尔卑斯复理石伊利石形态特征Sc=0.731，由这些关系式进行近变质带上下界线的转换得出与Kübler指数0.42～0.25。△2θ界线对应的Weaver指数和Weber指数的界线值分别为5.8～17.6和300～179。讨论了在自然干燥状态下直接进行伊蒙混层鉴定的基本原理。  相似文献   

基于TOUGHREACT的成矿过程化学反应数值模拟——以虎头崖铅锌多金属矿床为例          下载免费PDF全文

张婉秋  邹艳红 《地质找矿论丛》2020,35(3):354-362

虎头崖铅锌多金属矿床是青海省祁漫塔格成矿带中典型的矽卡岩型矿床。本文基于前人对该矿区地质特征、矿物学和地球化学特征等方面的研究成果，利用TOUGHREACT数值模拟软件，针对石英-晚期硫化物成矿阶段进行成矿过程化学反应数值模拟研究。通过建立地质概念模型和网格化物化参数的设置，在设置的时间步长内对多相流体流动和热对流传导、溶质运移以及地球化学反应过程进行模拟计算，研究了矿区内不同温压条件下成矿元素化学平衡浓度变化的总体趋势以及对方铅矿和闪锌矿矿物溶解沉淀情况的影响。结果表明：当温度从250℃降至155℃时，Pb²⁺的化学平衡浓度出现骤减；当温度从300℃降至100℃时，Zn²⁺的化学平衡浓度出现骤减；由于矿物的沉淀需要消耗大量相关的成矿元素，方铅矿和闪锌矿最佳成矿温度范围分别为155℃~250℃和100℃~300℃；不同温度和压力对成矿元素化学平衡浓度影响不同，温度是控制方铅矿和闪锌矿沉淀化学反应的核心因素，压力的影响较小。模拟结果与传统地质解译结果相对符合，证明了模拟的有效性。  相似文献   

橄榄岩中尖晶石的氧同位素继承性与晶体结构转变     

李一良  王峥嵘 《地球化学》1998,27(5):442-451

根据激光探针分析含尖晶石橄榄石中矿物的氧同位素组成,反应增量法计算矿物对氧同位素分馏曲线,得到正尖晶石－矿物对氧同位素温度高于含尖晶石橄榄岩的相平衡温度,而用反尖晶石－矿物的对计算的氧同位素温度则很接近含尖晶石橄榄岩的相平衡温度,较为合理,因此地幔橄榄岩中镁铝尖晶石的氧同位素组成继承了其母体反尖晶石特征,即与橄榄石,单斜辉石和斜方辉石平衡的是反结构尖晶石,在其发生相变作用变成正结构尖晶石时未发生氧  相似文献   

天山乌鲁木齐河源1号冰川物质平衡对气候变化的敏感性研究   总被引：42，自引：23，他引：19  

刘时银 《冰川冻土》1998,20(1):9-13

应用度日物质平衡模式对天山乌鲁木齐河源１号冰川物质平衡及平衡线高度对气候变化的敏感性进行了研究．结果表明,位于大陆性气候区且具有暖季补给特征的乌鲁木齐河源１号冰川物质平衡对气候变化的敏感性要小于海洋性冰川,升温１℃或增加２０％的降水可引起平衡线上升８１ｍ或下降３１ｍ．此外,气温与降水在物质平衡形成过程中的作用是不同的,气温引起物质平衡剖面以旋转方式变化,而降水可导致其平移方式的响应．若未来升温２℃时,即使降水增加３０％,１号冰川向负平衡变化仍然不能得到遏制．  相似文献   

Theoretical investigation of large-amplitude waves in granular soils     

V.A. Osinov 《Soil Dynamics and Earthquake Engineering》1998,17(1):13-28

A theoretical investigation of plane waves in granular soils is presented. Dynamic equations are derived with the use of the hypoplasticity theory for granular materials. For numerical calculations the material parameters of Karlsruhe sand are used. Wave speeds as slopes of characteristics of the dynamic equations are calculated for various stresses and densities. It is shown that under certain conditions the dynamic equations lose hyperbolicity and the initial boundary value problem thus becomes ill-posed. Two types of ill-posedness are found, known as flutter ill-posedness and stationary discontinuity. The latter is shown to arise at higher shear stress than the former. A comparison is made between dynamic ill-posedness and stability of static equilibrium. With the use of the second-order work stability criterion it is found that the dynamic equations lose hyperbolicity when the static equilibrium under a dead load is still stable. Numerical solutions to the problem of propagation of boundary disturbance in a half-space are obtained. Owing to dilatancy and contractancy of the granular material, a purely transverse disturbance induces a longitudinal component of velocity in the wave, and vice versa.  相似文献   

藕塘古滑体在三峡水库形成后的稳定性分析   总被引：2，自引：0，他引：2  

杨家岭  朱维申  罗晓东  张玉军 《岩土力学》1998,(2)

采用极限平衡和三维弹塑性有限元分析方法,对藕塘古滑体在三峡水库形成后的稳定性进行了计算分析。计算中考虑了水库水位变化、渗透压力、建筑荷载以及地震等诸多因素对滑体稳定性的影响,并给出了滑体失稳的临界破坏条件,为库区移民迁建新址的勘察工作提供了量化的论证。  相似文献   

Nonlinear Coupled Analysis of a Single Point Mooring System     

WANG Junrong LI Huajun LI Ping ZHOU Kai 《中国海洋大学学报(英文版)》2007,6(3):310-314

Coupled effects on a single point mooring(SPM) system subjected to the combined action of wind,waves and current are studied in this paper. Due to the complicatedness of the sea state and the huge size of the vessel,physical experimental study is both time consuming and uneconomical,whereas the numerical study is cost-effective and DNV software provides powerful SESAM software in solving the issues. This paper focuses on the modeling process of the SPM system,catenary equilibrium calculation,static analysis of the vessel in three different scenarios,and dynamic response simulation of the SPM system under environmental excitations. The three scenarios in study are as follows:the SPM is under the combined function of(a) wind,waves and current,(b) wind and waves,(c) current and waves. They are so set that one can compare the contributions of different types of loads in both static and dynamic studies. Numerical study shows that wind and current are the two major factors contributing to the mooring line tension,and surge and sway are the two dominant motions of the moored vessel subjected to environmental excitations.  相似文献   

An enlarged and updated internally consistent thermodynamic dataset with uncertainties and correlations: the system K2O–Na2O–CaO–MgO–MnO–FeO–Fe2O3–Al2O3–TiO2–SiO2–C–H2–O2     

T. J. B. HOLLAND  R. POWELL 《Journal of Metamorphic Geology》1990,8(1):89-124

We present, as a progress report, a revised and much enlarged version of the thermodynamic dataset given earlier (Holland & Powell, 1985). This new set includes data for 123 mineral and fluid end-members made consistent with over 200 P–T–X_CO2–f_O2 phase equilibrium experiments. Several improvements and advances have been made, in addition to the increased coverage of mineral phases: the data are now presented in three groups ranked according to reliability; a large number of iron-bearing phases has been included through experimental and, in some cases, natural Fe:Mg partitioning data; H₂O and CO₂ contents of cordierites are accounted for with the solution model of Kurepin (1985); simple Landau theory is used to model lambda anomalies in heat capacity and the Al/Si order–disorder behaviour in some silicates, and Tschermak-substituted end-members have been derived for iron and magnesium end-members of chlorite, talc, muscovite, biotite, pyroxene and amphibole. For the subset of data which overlap those of Berman (1988), it is encouraging to find both (1) very substantial agreement between the two sets of thermodynamic data and (2) that the two sets reproduce the phase equilibrium experimental brackets to a very similar degree of accuracy. The main differences in the two datasets involve size (123 as compared to 67 end-members), the methods used in data reduction (least squares as compared to linear programming), and the provision for estimation of uncertainties with this dataset. For calculations on mineral assemblages in rocks, we aim to maximize the information available from the dataset, by combining the equilibria from all the reactions which can be written between the end-members in the minerals. For phase diagram calculations, we calculate the compositions of complex solid solutions (together with P and T) involved in invariant, univariant and divariant assemblages. Moreover we strongly believe in attempting to assess the probable uncertainties in calculated equilibria and hence provide a framework for performing simple error propagation in all calculations in thermocalc, the computer program we offer for an effective use of the dataset and the calculation methods we advocate.  相似文献   

P–T evolution of kyanite eclogite from the Pirin Mountains (SW Bulgaria): implications for the Rhodope UHP Metamorphic Complex     

M. JANÁK  N. FROITZHEIM  N. GEORGIEV  T. J. NAGEL  S. SAROV 《Journal of Metamorphic Geology》2011,29(3):317-332

A new occurrence of kyanite eclogite in the Pirin Mountains of southwestern Bulgaria within the rocks belonging to the Obidim Unit of the Rhodope Metamorphic Complex is presented. This eclogite provides important information about the peak–pressure conditions despite strong thermal overprint at low pressure. Textural relationships, phase equilibrium modelling and conventional geothermobarometry were used to constrain the metamorphic evolution. Garnet porphyroblasts with inclusions of omphacite (up to 43 mol.% Jd), phengite (up to 3.5 Si p.f.u.), kyanite, polycrystalline quartz, pargasitic amphibole, zoisite and rutile in the Mg‐rich cores (X_Mg = 0.44–0.46) record a prograde increase in P–T conditions from ～2.5 GPa and 650 °C to ～3 GPa and 700–750 °C. Maximum pressure values fall within the stability field of coesite. During exhumation, the peak–pressure assemblage garnet + omphacite + phengite + kyanite was variably overprinted by a lower pressure one forming symplectitic textures, such as diopside + plagioclase after omphacite and biotite + plagioclase after phengite. The development of spinel (X_Mg = 0.4–0.45) + corundum + anorthite assemblage in the kyanite‐bearing domains at ～1.1 GPa and 800–850 °C suggests a thermal overprint in the high‐pressure granulite facies stability field. This thermal event was followed by cooling at ～0.8 GPa under amphibolite facies conditions; retrograde kelyphite texture involving plagioclase and amphibole was developed around garnet. Our results add to the already existing evidence for ultra high pressure (UHP) metamorphism in the Upper Allochthon of the Rhodope Metamorphic Complex as in the Kimi Unit and show that it is more widespread than previously known. Published age data and field structural relations suggest that the Obidim Unit represents Variscan continental crust involved into the Alpine nappe edifice of the Rhodopes and that eclogite facies metamorphism was Palaeozoic, in contrast to the Kimi Unit where age determinations suggest a Jurassic or Cretaceous age for UHP metamorphism. This implies that UHP metamorphism in the Upper Allochthon of the Rhodopes may have occurred twice, during Alpine and pre‐Alpine orogenic events, and that two independent HP/UHP provinces of different age overlap in this area.  相似文献